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Abstract 

We consider the equation that models the spreading of thin liquid films 
of power-law rheology. In particular, we analyze the existence and unique- 
ness of source- type self-similar solutions in planar and circular symmetries. 
We find that for shear-thinning fluids there exist a family of such solu- 
tions representing both finite and zero contact angle drops and that the 
solutions with zero contact angle are unique. We also prove the existence 
of traveling waves in one space dimension and classify them. 



1 Introduction 

Here we study capillary spreadings of thin films of liquids of power-law rheology, 
also known as Ostwald-de Waele fluids [8] . The power-law rheology is one of the 
simplest generalizations of the Newtonian one, in which the effective viscosity rj 
at a point is assumed to be a function of the local rate of deformation 7 given 
by t] = ml'yl 1 /* -1 . The values of m and A depend on the physical properties of 
the liquid. When A > 1 the viscosity tends to zero at high strain rates [8] and 
is larger at low strain rates (these fluids are called shear-thinning). 

In [10, 6] the following equation for one dimensional motion was derived 
using the lubrication approximation: 

u t + (u x+2 \u xxx \ x ~ 1 u xxx ) x = 0. (1) 

where u(x,t) represents the thickness of the one-dimensional liquid film at po- 
sition x and time t. In [1] a generalized version of (1) was studied by means 
of asymptotic and perturbative techniques in order to construct approximate 
solutions representing the spreading of a droplet. We will look for solutions 
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with a compact support [x\(t), X2(t)]. Therefore, by conservation of mass and 
equation (1), 

= — / u(x, t)dx 
dt J Xl (t) 

= x' 2 (t)u(x 2 (t), t) - x[(t)u(xi(t), t) + (u^^xxx^^Uxxx) 

which forces us to impose 

x' l {t)u{ Xl {t),t) + (u A+2 |Usxx| A_1 U*s*) { Xi (t),t) = , i = 1,2. (2) 

Equation (1) formally admits solutions of the form 

u{x,t) = ^u{^) (3) 

with 

^=5AT2- < 4 > 

Without loss of generality, one can set A 2X+1 = l/(5A+2). Then with tj = x/t 13 , 
U (rj) satisfies the ODE, 

U X+2 U" ,X = V U , (5) 

which results from introducing (3) into (1), integrating once and choosing the 
integration constant K — (by(2)). 

The initial conditions for symmetric drops are 

U(0) = 1, U'(0) = 0, U"(0) = -k, (6) 

where k is a real positive parameter. 

For radially symmetric flows, the equation of motion is (cf. [7]): 



u, 



+ - [ru x+2 (u rr + — ) (u rr + —)) = 0, 



where r is the radial coordinate, and the self-similar solutions are of the form 

J)l f 7* \ 

<r,t) = [^) 

with 

1 



7A + 3 

and U(r]) satisfying the following ordinary differential equation: 

U x+2 Uu" + jj) =vU. (7) 
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Analogously to (6) we impose 

1/(0) = 1, U'(0) = 0, lim -(nU')' = -K. (8) 

For simplicity and clarity we begin with a datached analysis of the existence 
and uniqueness of one-dimensional self-similar solutions, i.e., solutions to Eqs. 
(5-6). This is done in section 2. In section 3 we indicate the changes needed 
to deal with radially symmetric self-similar solutions (Eqs. (7-8)). Finally, in 
an Appendix we study the traveling wave solutions of (1) and classify all the 
possible behaviors of moving fronts close to the interface. 



2 Analysis of one dimensional self-similar solu- 
tions 

In order to remove the parameter k from the initial condition we introduce the 
change of variables 

x — tj^/k. (9) 

Then z(x,j) = U(rj) satisfies 

z 1+a z'" = 1X a (10) 
z(0) = l,z'(0) = 0,z"(0) = -1, 

with 

7 = fc-(3+l/A)/2 (n) 

and 

o= 1/A < 1. 
In the following subsections we prove: 

Theorem 1 For each a € (0, 1) there exists a 7 = 7(a) and y = y(a) such 
that z(x,-f(a)) first reaches z = for x = y(a) and z'(y(a), 7(a)) = 0. The 
function z(x,j(aj) satisfying these conditions is unique. 

Theorem 2 Given a € (0, 1) and 8 S (0, 1) there exist 7 = 79 and y = ye 
such that z(x, 70) first reaches z = for x = ye and z'(ye,"/g) = —y/28. 

These results are physically meaningful, because they imply that for A > 1 
there exist solutions describing one-dimensional drops with fronts advancing at 
finite speed. This result contrasts with the Newtonian case A = 1, where such 
solutions do not exist [3]. 

We also show that source- type self-similar solutions cannot exist for A < 1. 

All these results were suggested by numerical calculations in [6] and asymp- 
totic analysis in [6, 10]. 

In an appendix we will show, by studying traveling wave solutions, that the 
local behaviors of moving fronts near the interface are just those described in 
Theorems 1 and 2. 
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2.1 General Properties of Solutions. 



We consider the initial value problem 

z l+a z m = ^ 
z(0) = l,z'(0) = 0,z"(0) = -1, 

where a > is fixed and 7 £ R is a parameter. We are interested in finding those 
values of 7 for which there are interfaces, i.e., such that the solution z = z(x, 7) 
has a zero for a finite value of x and 7. 

If 7 = then the non-trivial solution to (12) is 

*(s,0) = l-y (13) 

which decreases from 1 to as a; increases from to The interface is x = 

If 7 < then the first three ^-derivatives of z(x, 7) are negative whenever z > 
0. Thus z is a decreasing function of a: as long as z is positive. Moreover, z'" < 
and the initial conditions imply that z(x,j) < z(x,0) for x > 0. Physically, the 
case 7 < represents retracting, dewetting drops. 

If 7 > then z(x, 7) > z(x, 0) for x > 0. Moreover, z'(x, 7) < at least for 
all sufficiently small x > 0. Either z'(x,j) < for all x such that z(x,j) > 
or else there exists a ( > such that z(x,j) decreases on (0,() an d achieves 
a minimum at x = ( with z(x,j) e (0, 1). At x — ( we have z'((,-f) — and 
z"((,j) > 0. Since z'" > it follows that z"{x,i) > and z'{x,i) > for all 
x > (. Therefore if z is not everywhere decreasing for z > then z has a unique 
minimum at some x — ( > with 

< z(C, 7) < z(x, 7) < 1, z'(x, 7) < for x <G (0, C) 

and 

2(C,7) < z{x,l), z'(x,i) > for xe (C,oo). 



2.2 Local Expansions Near the Interface 

There are two different asymptotic expansions for the solution near the root 
z(y,-f) = 0: one represents solutions with finite contact angle, 

z(x, 7) = V29(y - x) H j— — (y — x) 2 ~ a + lower order terms 

2— e 1+a a{\ - o)(2 - a) 

(14) 

where G (0, 1), and another with zero contact angle 

/ icy , V3 \ V(2+0) 

z ( x?7 ) = + a ^ ] y a /( 2 + a )(y-x) 3 /( 2+Q )+lower order terms. 

\3a(l — a)(l + 2a) / 

(15) 
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It is simple to show the nonexistence of local expansions near the interface 
if a > 1 and 7 > 0. If a > 1 and 7 < 0, then the local expansion is given by 
(14). 

We can write formal representation formulae for the solution to (12). Let 



zi+°(£,7) 



h (x, 7) = 
and 

i j+i( x ,i) = / ^(CiMC 

Jo 

for j = 1, 2. Then as long as z > 

z"(x,7) = -l + 7 / 1 (x,7), (16) 
z'(a;,7) = -a; + 7/2(2;, 7), (17) 

and 

x 2 

z{x,i) = l- — +'yh{x,i). (18) 

It is clear from (14) and (15) that if a < 1 then the representation formulae 
(17) and (18) can be extended to the interface. In particular, if z(x, 7) > for 
< x < y and z(x, 7) — > as x y then 

z"(x, 7 )^+oo, (19) 

z'(x, 7) -» z'(y, 7) = —j/ + 7/ 2 (y, 7), (20) 

and 

V 2 

z{x,i) -> z(y,7) = 1 - y +7/3(^,7). (21) 

If a > 1 then (19) and (21) are still valid, but (20) is replaced by z'(x, 7) — ► +00. 
Hence, there are no interfaces for a > 1 unless 7 < 0. 

There are no interfaces for sufficiently large values of 7. More precisely, we 
prove 

Lemma 1. There are no solutions with an interface if 

l + q 

Proof. As we observed above, any positive local minimum of z(x, 7) is in 
fact a unique global minimum. Thus if z has a positive minimum, then z docs 
not have an interface. Suppose z has a positive minimum at x = xq. Then since 
< z < 1, we have 

^ = 1 - T +7/3( *' 7) * 1 - T + (l + a )( 7 2Ta)(3 + fl) S P( ^ 7) 

for x € [0, xo]. The condition (22) guarantees that the positivity of the minimum 
value of P. Therefore if (22) is satisfied, z is everywhere positive and has no 
interface. □ 
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2.3 Analysis for Small 7 

Since (12) is singular at z = 0, we shall first study the behavior of the solutions 
for z > 0, and, instead of studying the interfaces z(y,-y) — directly, we first 
solve 

z(y,7) = 5 (23) 

for given 6 e (0, 1). In the next subsection we study the limit as 6 \ 0. We 
begin by studying (a) in the neighborhood of the solution (13) for small |7|. 

We have 

z{yi(S),0) = 6 and z'{ yi {S),0) = -yi(6) < 0, 

where yi(5) = y/2(l-5). Since z(x,j) is smooth for z > 0, by the Implicit 
Function Theorem there exists a smooth function 2/(7, S) for sufficiently small 
I7I such that 

V(0,S) = yi(6) and z(y(j, S), 7) = 6. 
2/(7, S) exists for all 7 such that 

*'(»(7,*),7)<0. (24) 

As we observed above, (24) holds for all 7 < 0. Moreover, since (24) holds for 
7 = 0, and since 2/(7, 5) depends smoothly on 7 and z (y, 7) depends smoothly 
on y it follows that (24) continues to hold for all sufficiently small 7 > 0. 
Differentiating (23) with respect to 7 we find 



z (y>Vg^ + z t(y> 7) = 0. 



Therefore 



dy _ z 7 M<5),0) J 2 (yiffl,0) ^ n 
9 7 ' 7 -° (5), 0) ^2(1^) 

so that j/ is an increasing function of 7 at least for I7I small. 



2.4 Estimates for solutions with interfaces 



We now investigate the maximal 7- interval of existence for the functions 2/(7, 5) 
constructed above. In the limit as S — > this will lead to solutions with zero 
contact angle. 

In view of (24) we define 

7o(J) = sup {7 > : 2' (2/(7, 8), 7) < 0} 

and 

2/0 (<5) = lim 2/(7,5). 

7/7o (<5) 
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Note that since z'{y,~f) and 2/(7, 5) are continuous we have 

z(y (5),jo(S)) = 5 and z'(y (S), l0 (S)) = 0. 
We now derive some estimates for 70 (5) and yo (<*>)• 
Lemma 2. Let 

fl(7) = { (1 + Q) 7 (2 + a) p and G = 2^ (1 + a)(2 + a). 

Then 

< 7o (5) < G (25) 

and 

v/2(l - 5) < yoW < S( 70 (5)). (26) 

Proof. For < x < 1/(7, 5) we have 5 < 2(2,7) < 1> an d z '{y{li *)> < 0- 
It follows from (20) that 

> «'(y(7,*),7) = -1/(7,*)+ 7/2(1/(7, «), 7) > -2/(7,*) + (/f^+i) ■ 
Therefore 

^2(1 - 5) < y(7,5) < 5(7). 

Letting 7 /" 70 (5) we immediately obtain (26). Then (25) follows from ^/2(1 — 5) < 
B( l0 (5)).D 

For arbitrary fixed 6 € (0, 1) consider 



A( 7 , 9) = z'(y( 7 , 5), 7) + 9^2(1-5). 

A(7, 0) is a continuous function of 7 e [0, 7o(<5)] with 

A(0, 9) = -(1 - 9)^2(1-5) < and A( 7o (5), 9) = 9^/2(1-5) > 0. 

By the Intermediate Value Property of continuous functions there exists 7(5, 9) e 
(0, 70 (5)) such that A( 7 (5, 9), 9) = 0, i.e., 

z'(Y(5, 6), 7(5, ^)) = -0 v /2(l — 5), (27) 

where 

r(M) = i/(7(<*, *),<*)• 

We take 7(5,0) to be the smallest positive value of 7 for which (27) holds. We 
extend the definition of 7(5,0) by setting 7(5,0) = 7o(<5) and 7(5,1) = 0. In 
view of (25) 

7 = lim inf 7(5, 9) 

<5\0 

exists and satisfies 

< 70 < G 
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for arbitrary £ [0, 1]. For each £ [0, 1] fix a sequence {5g} such that S J 9 \ 
and 7^ = ^{8 3 , 9) -> 70 as j -» oo. 

Lemma 3. 79 > if and only if a < 1. 

Proof. From (17) we have 

-9^2{l~S) = z'(Y(5, 9), 7 (5, 0)) - -F(5, 0) + 7(5, 0)7(5, 0) 

or 

F(5, 0) - V / 2(T~ 7 5) = 7 (5, 0)7(5, 0), 

where 

7(5,0) = I 2 {Y (6,9), 7(5,0)). 
Suppose that jg = 0. Then since 

(1- 0)^2(1 - 51) < ^1(51,9), 

7(5g, 0) — ► implies 7(5^, 0) — ► 00. In view of the asymptotic behavior of z near 
this can only occur if a > 1. On the other hand, if 7(5g, 0) > then 

7(5,0)7(5,0) <Y(S,9) <7?( 7 (5,0)) 

so that 
Therefore 

< lim sup 1(61,0) < < 00 

j^oc 70 

which can occur only if a < l.D 

We now investigate the behavior of Y(Sg, 0) as j — > 00. 
Lemma 4. If a < 1 £/ien 

V2<y e = lim K(<5*, 0) < B{ le ). 

j— >oo 

Proof. 7^ — > 7e implies that for arbitrary e £ (0,70) there exists an N — 
N(e) > such that 7^ > 79 - e for all j > N. In view of (26) 

^2(1 -Sl)<Y(Sl,9)<B( le -s). 

Thus there is a subsequence of the 8 3 e which we again call {8g} along which 
Y(5g, 6) converges to a limit ye which satisfies 

V2 < y 9 < B{ le - e). 

Since e is arbitrary the assertion now follows by letting e \ 0.D 
Finally we prove 
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Lemma 5. If a < 1 then 
z{y e ne) = Km z(F(^,0), 7e ) = and z'(y e , le ) = Hm z'(y(^, 0), 7e ) - -6s/2. 

J — *oo J — *oo 

Proof. From the representation formulas (17) we have 

z(Y(5i, 6), 7fl ) = 1 - Y2( f &) + loh(Y(Si,9), le ) 

which we rewrite in the form 

z{Y{5i,6), le ) = Si + le {h{Y{5 j e ,e), l0 ) - I 3 (H<&<M)} 
+ {le-l i e}h{Y{8le), 1 i). 

The first and third terms on the right hand side clearly tend to zero as j ; — > oo. 
The convergence of the middle term to zero follows from the asymptotic behavior 
of z near zero and the Lebesgue Dominated Convergence Theorem. The result 
for z' is proved in a similar manner. □ 
Corollary. If a < 1 then y g > 

Proof. Since -fg > we have z"'(x,jg) > on (0,ye). Therefore 

O = z{y 0lle ) > 1- y 
which implies that y$ > >/2.D 

The analytic considerations given above establish the existence for a < 1 
of solutions whose interfaces have zero contact angle and of solutions whose 
interfaces have non-zero contact angles. In particular, for each 9 £ [0, 1] there 
is a pair (ye,je) such that the solution z(x, r ye) has its interface at x — ye with 
z'(ye,7e) = —0V%- However, the results so far do not give any global informa- 
tion about uniqueness and monotonicity of (ye, 70). In the next subsection we 
will prove the uniqueness of zero contact angle solutions for each 7 > 0. 

2.5 Uniqueness of zero contact angle solutions 

In this section we present a proof of uniqueness of the zero contact angle solution 
following arguments similar to those used in [4] in a somewhat different context. 
Rewrite equation (5) in the form 

U'" =nilJ- s , (28) 
where s = 1 + 1/A, and let Uj{rf) for j = 1,2 be solutions to (28) such that 

^•(0) = 1,^(0)=0, 
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and 

Without loss of generality we assume that 
Set 

3A + 1 

Vi\ 2X+1 , ,_ Wi 



c = — and £ = 

Note that c > 1. The rescaled function 

is a solution to (11) with 

t7 2 (o) = c,^(o) = o, 

and 

u 2 (vi) = u' 2 {^) = 0. 

For £ G [0, ?7i] define 

K(0 = c/ 2 (0 - c/i(C). 

Then 

F(0) = c - 1 > 0, y' (0) = 0, and 1/(77!) = 0. 
Moreover, in view of (28), 

V"' = $ {/(%)- /(C/i)}, 

where 

/(t) = t" s . 

Note that / is a monotone decreasing function with /' < for r > 0. By the 
Mean Value Theorem 

W" = £*f'{W)V 2 
for some W between U2 and U\. Thus 

W" < on [0,771). (29) 

Set 
Then 

h(0) = (c - 1)V"(0) and ti = W" < on [0, 771), 

i.e., ft, is a monotone decreasing function on [0,771]. Using the asymptotic ex- 
pansion (15) we find that for sufficiently large £ < 771 
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for some constant K. Since a < 1 it follows that 

h(Vi) = 0. 

If h(0) < then h decreasing and /i(?7f) = imply that ft e on [0, 771] . 
Thus 

h' = ^(Uz 3 - Ui S )(U 2 -Ui)=0 on [0,771) 

and this implies that U2 = U\ on [0, 771). In particular, c = 1/2(0) = U\(0) = 1. 
It follows that iji — i]2 so that 

U 2 (0 = U2(0 = Ui(0 on [0,7?!]. 

Suppose that h(0) > 0. Together with h(rji) = and the fact that h is 
decreasing this implies that h > on [0, 771]. Therefore 

VV" > ^V 12 > on [0,771]. (30) 

Moreover, it follows from h(0) > that c > 1 and V"(0) > 0. Hence there exists 
a(i£ (0,771) such that 

V > c - 1 > on [0, Ci] and V"(Ci) > 0. 

On the other hand, V(r)i) = means that V, which is initially increasing, must 
start to decrease somewhere in (£1,771). Thus there exists a ( 2 € (Cb?7i) sucn 
that 

V > on [0,C 2 ] and V"(C 2 ) < 0. 
It follows that there exists a C3 € (Ci, C2) such that 

V(( 3 )V"(( 3 ) < 

and this contradicts (30). We conclude that we must have h(0) < and, as we 
have shown above, this implies that U\ = JJ^- 



3 The radially symmetric case 

In this section we will briefly describe how the proof of existence and uniqueness 
for the 1-D model extends to radially symmetric solutions. In this case, the 
evolution equation is (cf. [7]): 

/ 11 \ A-l 

A+2 



u t + - \ ru 

r 



and the self-similar solutions are of the form 

Al f t 
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with 

/3= 7X + 3 

and U(r]) satisfying the following ordinary differential equation: 

U x + 2 ((u" + — X) =r,U . (31) 



V 



Analogously to (6) we impose 



1/(0) = 1, U'(0) = 0, lim -(r]U')' = -K. (32) 

After a suitable rescaling, the problem (31), (32) becomes 

z 1+a (z" + -) = ix a , (33) 



z(0) = l,z'(0) = 0, lim -(ajz')' = -l- (34) 
x^o+ x 

We will denote the solutions to this problem by z(x,j). One can easily check 
that the asymptotics near the contact line are identical to the one dimensional 
case, and in particular, for the zero contact angle solutions are given by Eq. 
(15) (cf. [7]). When 7 = 0, there exists an explicit solution to (33), (34) given 
by 

x 2 

z(x,0) = l--. 

If 7 < then it is simple to show that z', z", z'" are all negative whenever 
z > 0, implying the existence on compactly supported solutions to (33), (34) in 
this case. 

If 7 > then one can write 

1 



(xz'y 
x 



A+a ■ 



Noticing that 



dz 

a; v ' ) v a dy 2 \^ dy / 

where y — x 2 it is clear that function yjj^ is concave up. At the origin yj^ = 

and ^ (j/^f ) < 0- Hence, ^ has at most one zero for y > 0. This implies 
z(x,-f) can either decrease monotonically to zero or grow monotonically after a 
minimum. 

The use of implicit function theorem is analogous to the 1-D case and shows 
the existence of zero contact angle solutions as well as non-zero contact angle 
solutions. Since the proof follows the same lines, we omit the details. 
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The proof of uniqueness of the zero contact angle solution is similar to the 
1-D case except for the fact that the equivalent to inequality (29) is now 

V (1 {llV >)'\ ' = {V »> + Yl-YL )V <o. (35 ) 
\r/ ) T] r/ z 

Without loss of generality we can assume that V^O) > 0. One can easily 
verify that V"(0) = \ lim 7) ^ 0+ . 

Case 1: lim ?? ^o+ ' < 0. Using the asymptotic expansion (15) we can 
conclude that there exists an interval (772,771) where V(rj) and V"{rj) have a 
given sign and, moreover, V(rj)V{rj)" > in that interval. 

Let rf[ be a point of the interval sufficiently close to r\\. Let 77* be the 
inflection point closest to r)\. Given the regularity of V, since V"(0) < 0, 
V(0) > and V(r])V(r])" > for 77 e (772,771) such a point must exist and 
V"{7]*) = 0. 

We will show that V(rj*)V(r)i ) > for 77* < rf[ . The proof applies to 
the case in which V{vj{) > (the proof in the case V{rf[) < is analogous). 
Suppose that V(rj*) < 0. This implies the existence of some point 77** € (77*, 771) 
at which a local maximum of V is achieved and U" (77**) < 0. Hence there must 
exist and inflection point in (77**, 771) which contradicts the fact that 77* is the 
inflection point closest to 771. 

Since V" (i]*) is negative for 77 € (77* — e, 77*) and positive for 77 € (77*, 77* + s) 

we have that VV" is increasing at 77* and hence -£j(VV") — V(r]*)V"' (r]*) + 

V'(V*)V"(V*) = V{rf)V'"{r}*) > and V {ri*)V{ri*) < 0. This implies {V"'(r]*)+ 
Y ^2D- - ^P-Wir]*) > which contradicts (35). 

Case 2: lim^^ + ^2ZJ_ > 0. In the interval (0,771) there must exists a local 
maximum 772 for V and we can repeat the argument of the previous case looking 
for an inflection point at the left of 772 instead of 771. 



4 Conclusions 

In this paper we have established the existence of solutions representing the 
spreading of drops in a model for the capillary spreading of Non-Newtonian 
fluids of power-law rheology. The solutions depend on the rheology exponent A. 
Here we prove that when A > 1 then for a given mass of fluid, there exists only 
one solution with zero contact angle at the horizontal substrate and infinitely 
many solutions with a finite contact angle. Both the spreading rate and the 
height of the drop are power laws whose respective exponents depend on A. For 
the case when A < 1 such spreading solutions do not exist. Both results are 
valid for both planarly and radially symmetric drops. 

The results presented here pose a number of interesting questions: 

1. Does these solutions represent the intermediate asymptotics of the spread- 
ing of a compactly supported drop?. In other words, does the solution as 
t — > 00 depend on the initial shape of the drop?. 
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2. Is there a selection criterion between zero and non-zero contact angle 
solutions?. This is perhaps an stability problem that can be recast in 
the following way: which of these solutions is stable with respect to small 
perturbations of the initial data?. 



5 Appendix: the traveling wave solutions in 1-D 

In this section we study the existence and asymptotic behavior of traveling wave 
solutions. These solutions are relevant, since they allow the determination of 
the local behavior of moving fronts near the interface. They are of the form 

u(x,t) = f(x + ct) . (36) 

Substituting (36) in (1) we get the equation 

c/ + / A+2 |/^| A - 1 A« = ^ (37) 
where K is a real constant. The condition (2) forces us to choose K — 0. Since 

IA«I A_1 A« = -cf- 1 -" 

one has that > if c < and < if c > 0. Therefore, 

i 

and we can remove the factor |c| T by performing the change of variable / — > 

i 

|c| 2A+1 / • Given the translation and £ — > — £ reflection invariance of the equa- 
tion we can assume without loss of generality that the solutions representing 
fronts are defined for £ > and the front is located at £ = 0. We will also 
assume < and arrive to 

/«« = -r 1 -* • (38) 

Let us define 

* = / 

y = r a u 

with a = 4^ and (3 = together with the change of variable defined by 

jb-OJ+o)^ = (39) 

so that the solutions of (38) are orbits of the following third order autonomous 
dynamical system: 

x' = xy (40) 
y' = z - ay 2 (41) 
z' = -l + pyz (42) 
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'/ / / y / 
/ / / / / 
/ / / / / 
/ / / / / 
////// 



Figure 1: The phase plane for (41), (42) 



Notice that equations (41), (42) are uncoupled to equation (40) so that the 

problem reduces to the analysis of the phase plane for equations (41), (42). 

— - 1 2 \ 

There exists a unique equilibrium point at P : (y,z) — ((/3a) 3 ,a^d~^). 
This point is a saddle point as one can easily check. Hence, there exists an 
stable manifold Si through P as well as an unstable manifold £2. £1 consists 
on the union of two separatrices Ti, T2 and £2 on the union of another two 
separatrices T 3 , T 4 (see figure 1). There exist four different asymptotic behaviors 
for the trajectories as |y| tends to infinity and each of these asymptotic behaviors 
corresponds to the one exhibited by Ti, T 2 , T 3 , T4 respectively. Now we proceed 
to describe these behaviors in detail: 

1) Ti and T 4 are such that y — > ±00 and z — > 0* respectively. Then, from 
(41), (42) one obtains 



y ~ -ay , z ~ -1 + f3yz 



from which it follows 



dz 1 -fiyz _§_ 1 1 

— ~ = — z ~ Cy " — 7 — - ~ —7 -7— as y — > ±00 (43 

dy ay 2 y (a-p)y (a-0)y y V ; 

§_ _ A+2 
A-2 

2) r 3 and are such that z — ► +00 and y — » ±00 respectively. Then, from 
(41), (42) one obtains 

z' ~ 0yz , y' ~ z- ay 2 

from which it follows 

dz Byz , 3. , _2« /? 9 

^— ^ => (a + '-)y 2 ~dz a + z => z ~ (a + ^ j/ 2 as y -» ±00 (44) 

dy z — ay z 2 2 

since ^ > 0. 
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Once we have discussed the asymptotic properties of the trajectories in the 
phase plane corresponding to (41), (42), we proceed to study the behavior of the 
trajectories in the three-dimensional phase space for (40), (41), (42). Every man- 
ifold of the form KxT with V being a trajectory in the phase plane, is invariant. 
As we will see, the most interesting case corresponds to T = Tj (i = 1,2,3,4). 
The invariant manifolds IT = M x I\ intersect at the points IxP which form 

the trajectory (x,y,z) — (e^ Q ) 3 ^ 1 , (/3a)^^ ,ai/3~%). The behavior of the tra- 
jectories on IT (i = 1,2,3,4) arc rather different. The trajectories in IT and 
±I 2 approach asymptotically RxPasn +oo. From equations (40), (41) and 
the asymptotic behavior computed in formula (43) it follows that y ~ Kx~ a 
as x — > + (K an arbitrary positive constant) for the trajectories in IT. Anal- 
ogously, from equations (40), (41) and the asymptotic behavior computed in 
formula (44) it follows that y ~ Kx? when y — > — oo (K an arbitrary negative 

constant) for the trajectories in IT. The trajectories in II3 and II4 start at 
_ i 1 2 

(x,y,z) = (0, {f3a) 3 ,a3p~i). Analogously to the trajectories in n 2 and IT, 
the trajectories in 1I3 are such that y ~ Kx^ when x — > 00 (K an arbitrary 
negative constant) and the trajectories in II4 are such that y ~ Kx~ a as x — > + 
(K an arbitrary positive constant). 

Finally, we translate the phase space trajectories described above into solu- 
tions of (38) and discuss their physical significance. The trajectory in M x P is 
such that, by (39) and (40), we have 

^=x 1 ~^ (/?<*)"* 

which implies, imposing x(0) = 0, that 

/(0 - x(0 = [(a + 0) (f3ayi] ™ ^ = C A £*£r . 

The trajectories in IT are such that /(£) <~ CaC 2A+1 as C — ► 00 and, given 
that y ~ Kx~ a as x — > + , by (39) and (40) one has 

so that /(C) ~ K£ as C ^ 0+. 

The trajectories in II 2 are such that /(£) ~ CaC 17 ^ 1 as C — > +00 and 
/(C) ~ K£ 2 as C — > —00. Since there is no front in this case, we do not consider 
this solution to be physically relevant. 

The trajectories in II3 are such that /(C) ~ CaC 2 -*^ 1 as £ — > + and /(C) ~ 
i^C 2 as C — > 00. 

The trajectories in II4 are such that /(C) ~ CaC 2A+1 as C — » + and /(C) ~ 
if (Co — C) as C — » Co" f° r some positive Co- They are compactly supported. 

In addition there exist another four families of solutions such that the tra- 
jectories in phase space approach asymptotically to two of the manifolds IT 
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(i = 1,2,3,4). The first one approaches 111 and H4 so that behaves linearly 
at the origin and linearly close to some £0 > 0. It is compactly supported and 
lacks a clear physical significance. The second one approaches II4 and U 2 so 
that it is linear at the origin and grows quadratically at infinity. For them y < 
which implies that /'(£) < 0. They represent dewetting solutions. The third 
one approaches n 2 and II3 and grows quadratically at ±00 presenting no fronts. 
The last one approaches II3 and IT. It behaves linearly at the origin and grows 
quadratically at infinity. 

To summarize, one has infinitely many solutions that behave linearly at the 
origin and such that /(£) ~ Ca£ 2A+1 as £ — > 00. There is only one solution 
(which is explicit) such that /(£) ~ Ca£ 2A+1 at the origin (zero contact angle) 
and at infinity. There exist infinitely many solutions with zero contact angle at 
the origin and growing quadratically at infinity and, finally there exist dewet- 
ting solutions which are linear at the origin and grow quadratically at infinity. 
Hence, the only local behaviors near the interface for moving fronts are those 
with finite contact angle and the one with zero contact angle. 
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